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Abstract 



> 

£^ ' In this paper the second Lyapunov method is used to study the 

, stability of the de Sitter phase of cosmic expansion when the source of 

c"^i ' the gravitational field is a viscous fluid. Different inflationary scenarios 

0^ , related with reheating and decay of mini-blackholes into radiation are 

investigated using an effective fluid described by time-varying thermo- 

O dynamical quantities. 

a -1 

60- 1 Introduction 



In recent years considerable attention has been paid to the bulk-viscosity 
driven inflationary scenario. This is only natural since the effect of bulk 
viscosity in an expanding universe is to reduce the equilibrium pressure. 
Therefore one may wish to know whether this effect could be strong enough 
to render a large negative effective pressure that leads to inflation. 

Fundamental strings can create an initial cosmological state either of 
exponential or power law inflation followed by a smooth evolution towards 
the typical Friedmann decelerated expansion dominated by the stress-energy 
tensor of a radiation fluid. This connection is natural because the string- 
driven inflationary expansion may arise due to the spontaneous quantum 
production of fundamental strings on scales larger than the horizon, and 
a phenomenological bulk viscosity can be used to describe the effect of par- 
ticle production 0,0. For some cosmological implications of fundamental 
strings see also [1 and references therein. 
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Inflationary scenarios have usually been associated with the dynamics of 
a spatially homogeneous scalar field such that its potential energy overpow- 
ered the kinetic energy, and the equation of state the of vacuum p^ = —p^ 
was satisfied. If at the time of interest the scalar field dominated any other 
form of energy, then the cosmic scale factor increased exponentially with 
time ||. The fact that a(t) > 0, where a(t) is the scale factor, has several 
nice consequences, among others a scale invariant spectrum of initial density 
perturbations, the one more likely to be compatible with observation |J. 

However inflation, either exponential or power-law, can in principle be 
driven by any mechanism that renders the total hydrostatic pressure nega- 
tive, such as bulk viscous pressure associated with non-adiabatic expansion 
in FLRW universes -the effectiveness of this mechanism has been discussed 
in the literature @. Usually in these models a mixture of relativistic and 
nonrelativistic (heavy) particles is assumed at some early phase of cosmic ex- 
pansion. There the bulk viscosity can be very large and may drive inflation. 
After the decay of these heavy particles bulk viscosity vanishes, terminating 
the inflationary phase, and returning to the radiation dominated Friedmann 
universe ||. 

Very often a viscous pressure represents only a small perturbation to the 
equilibrium (hydrostatic) pressure of the fluid. However, as is well known 
-see H and (T^j- the effect of particle decay can be phenomenologically un- 
derstood as a dissipative pressure, and this one can be very large depending 
on how big the decay rate is. This approach was put on a solid footing by 
Triginer et al. 11]. 

Stability of the de Sitter solution in cosmological models where the source 
of the metric is a dissipative fluid obeying some causal transport equation 
has been investigated by several authors. For the truncated version of the 
transport equation for viscous pressure, stability has been investigated in ||] 
and [p^| . In 13], |L4|], [15] and [ 16 1 the full version of the transport equation 
-see (||) below- was adopted. The nonlinear full version [jl7| was used in |18|| . 
The rationale behind causal transport equations can be found in 1£], for a 
short introduction see [20]. 

Here we introduce a new approach by using the second method of Lya- 
punov (2l]] to examine the asymptotic stability of the de Sitter solutions. 
This provides useful information on the dynamical behavior of the system, 
not only near the stationary solutions but far away from them as well. 

In section 2 the basic equations and Lyapunov's criterion for stability 
are presented and applied to different situations in which some or other of 
the three key parameters are held constant and the rest allowed to vary. In 
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section 3 the whole set of parameters are assumed to vary. Finally section 
4 summarizes the main findings of this work. 

2 Asymptotic stability 
2.1 General setting 

The spatially-flat Friedmann-Lemaitre-Robertson- Walker (FLRW) universe 
is described by the metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (1) 

Because of the spatial isotropy and homogeneity assumptions, velocity 
gradients causing shear viscosity and temperature gradients leading to heat 
transport are absent. The only possible dissipative term corresponds to the 
bulk viscosity (which as mentioned before may be interpreted as the effect 
of particle production). So the energy-stress tensor of this viscous fluid is 
given by 

T ab = pu a u b + (p + vr) h ab (2) 

where p is the total energy density of the cosmic fluid, p its equilibrium 
pressure, tt the dissipative scalar pressure, u a the four- velocity, normalized 
so that u a u a = — 1 and the tensor h ab = g ab + u a u b projects any tensorial 
quantity into the hypersurface orthogonal to u a . 
The corresponding Einstein field equations read 

3H 2 = np (3) 



H 



(p + p + vr) 



(4) 



where H = a/ a is the Hubble function, and k is Einstein's gravitational 
constant. We use units such that c = kg = h = 1, so that k = 8n/Mp, 
where Mp is the Planck mass. An over-dot denotes derivative with respect 
to cosmic time. 

As is well known, the dissipative pressure obeys the causal evolution 
equation |0] 



TT 



i + C 



2T( 



+ T7T 



-3(H 



(5) 
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where £ indicates the phenomenological coefficient of bulk viscosity, T the 
temperature of the cosmic fluid, and r the relaxation time associated to the 
dissipative pressure. Usually the latter is given by the kinetic theory of gases 



or by a fluctuation-dissipation theorem or both [24]. Provided the divergence 
on the left hand side is small the last equation can be approximated by 



7T + T7T 



(6) 



In view of equation (13) -see below- and the frequently made assumption 
that ( oc /), which is very natural for radiative-like fluids [23], pi], the 
above approximation may be justified for expanding regimes with H^^> H, 
something to be expected at the commencement of the Universe expansion. 
In any case the two above equations meet the requirements of causality and 
stability to be fulfilled by any physically acceptable transport equation ||2l 
At this point it is expedient to introduce the polytropic index 
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1 + P - 
P 



(7) 



which in general depends on time. This is obvious for a cosmic fluid consist- 
ing of a mixture of massless and massive particles, as these two components 
redshift at different rates and consequently 7 varies with time. 
This set of equations combines to 



-1 



H + 3jHH + t 



and the latter can be recast as 



H 



- (7 + r 7 ) H 2 



3 -CH 
2 s 







d_ 
~dt 



H 2 + V(H) 



D(H,H) 



(8) 



(9) 



Here, the left hand side is the time derivative of a Lyapunov function (see 
the appendix) with 



V(H) = t 



(10) 



and 



D(H, H) = - (3jH + r" 1 ) H 2 + \-t- 2 ( 



T7 



fry + r 2^) H 3 
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+ \r~ 2 (Cr - Cr) & (11) 

Assuming that the ratio £/r is bounded V(iif) has two extrema (see Fig. 
1), a maximum at = 0, and a minimum at 

#=^^- (12) 
7 + T7 



Inserting £ from ( |l2| ) in (|ll]) we get 



Q 

D = — (&yH + r" 1 ) if 2 - -t' 1 (7 + r 7 ) F 2 # 

- ^r" 2 (r7 - f 7 + r 2 7 ) # 3 (13) 

In the simplest case when the functions 7, r and Q are constants, D is 
seminegative definite in a neighborhood of the de Sitter solution H\ = 
C/ (7 + T i) an d the Lyapunov function has a time independent upper bound 
for large times, implying the stability of this solution On the other 

hand, when i7 is slowly varying we can neglect H in (|13|); thereby -D re- 
duces to 

^^"4^2 (r 7 - f 7 + r 2 7 ) # 3 (14) 

Provided D is negative definite, and taking into account that the Lyapunov 
function has an infinitesimal upper bound in the neighborhood of point 
(Hi,Q) of the phase space (H,H), the de Sitter solution fll2| ) is asymptoti- 



cally stable [21] 



In contrast with the usual inflationary scenarios driven by a large effec- 
tive cosmological constant necessarily accompanied by strong supercooling 
and reheating, dissipative processes in our scheme reheat the medium gently. 
In the following, several simple cases when either 7 or r remain constant are 
investigated. 



2.2 Variable r 

In the first instance (|l2] ) is asymptotically stable when r is a decreasing 
function, while in the second case stability occurs when 7 increases slowly. 
Thus, for instance, if 7 is a constant and r is a function that decreases 
in a first stage and then increases, we have that there is a first period of 
exponential inflation with H = C/7, and then a graceful exit. 
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For instance, using the relationship for the speed v of the dissipative 
signal 



v 2 jp (15) 



derived in [ 20 1 , we obtain 



(16) 



3v 2 H 

and therefore r will present a minimum provided v 2 has a maximum. In this 
case it may be said that the dissipative effect both drives an inflationary 
stage and then causes the exit from it. 

Another way to implement this scenario is as follows. Assume the cosmic 
fluid is modelled by a mixture of radiation and heavy particles (or massive 
modes of fundamental strings) that decay at a very high rate into (more 
stable) lighter particles (less massive modes), with high or moderate multi- 
plicity. Since the relaxation time for the interaction between radiation and 
massive particles is on very general grounds given by r ~ (no - ) -1 , where 
n denotes the number density of massive particles and a the interaction 
cross-section, which can be constant, there will be two stages. The first one 
corresponds to the decay of the more massive particles; there (despite the 
expansion) n will augment and r decrease accordingly. The second stage will 
commence when most of these particles have decayed; then the expansion 
will make n decrease and correspondingly r will increase ||25||. 



2.3 Variable 7 and decay of massive particles 

This second case will be useful to model the cosmic evolution during the 
inflationary period by allowing massive (dust) particles to decay into rel- 
ativistic particles [27]. Thus one may have simultaneously inflation and 



reheating. 

Exit from inflation occurs when particle production ceases and the fast 
dilution of relativistic particles makes 7 decrease again. An interesting 
model for the reheating-inflationary stage arises when r and £ are constant 
but 7 varies in such a way that Hi remains constant. Using ( |il"l) we see 
that D is semidefinite negative and therefore the de Sitter stage is stable. 
In this case 7(t) takes the form 

7 (t) = 7o + C 2 e-*/ T (17) 
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and H\ = C/7o- It is straightforward to implement a cosmological model in 
which 7 = 1 at t = t\ and 7 = 4/3 at t = £2, with £2 > h 



^) = 1 + 3(e -w T - e -.,/ T) < 18 > 

The evolution of 7(i) depends on the relationship between the timescales r 
and t2 — ti. When t <^ t 2 — t± the relativistic stage is reached very quickly. 
This corresponds to a universe initially dominated by very massive (dust) 
particles, that spontaneously decay into radiation at a high initial rate |27|. 
On the other hand, when r ^> t 2 — t\, the growth of 7(i) becomes milder. 

Since the detailed behavior of a dissipative relativistic fluid in a FLRW 
universe strongly depends on the thermodynamical properties of the cosmic 
fluid and generally these are poorly known, we are led to use the Boltzmann 
gas as the cosmological medium to obtain analytical results. Thus, while 
a Boltzmann gas is not a realistic model for the cosmological fluid in the 
actual Universe, it is a fluid for which the thermodynamical properties are 
well enough established by relativistic kinetic theory to allow us to build 
up precise models. The thermodynamic properties of the Boltzmann gas 
may be described by the dimensionless inverse temperature z = m/T, the 
relativistic chemical potential a and a constant Aq = m 4 / (2ir 2 ) , where m 
is the particle mass and g* the spin weight of the fluid particles. The ideal 
gas law has the form 

P = A e-^ (19) 
z z 



Ki(z) | z K 2 (z) 



z 2 



(20) 



where K n are modified Bessel functions of the second kind. Then, inserting 
(|l|) and © in ©, and assuming a vanishing chemical potential, we obtain 

l{z) = 1 + K ^ z ) ( 2i) 

R ' zK x {z) + 3K 2 (z) K ' 

Equating (|l8|) and (21) we can describe the continuous process of de- 
cay of mini-black holes from t = t\ when the black hole energy density 
dominates the Universe, until t 2 when the black holes have completely evap- 
orated away and the Universe is radiation-dominated. In so doing we are 
implicitly assuming that all the black holes have the same mass and there- 
fore the same temperature, and that this one equals the temperature of the 
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massless component of the cosmic fluid at the beginning of the evaporation. 
Scenarios in which this situation may occur have been reported in the lit- 
erature -see for instance ||2S|| , [ pS| ] and [ j3C|| . In |2^] black holes are formed 
as consequence of the first order phase transition from the false to the true 
quantum vacuum; in |2^] because of collisions between bubbles of the new 
phase, and in |3(J by quantum fluctuations at the end of hybrid inflation. 
In all three cases these abundantly produced mini-black holes dominate the 
Universe and their subsequent explosive evaporation into lighter particles 
can be modeled ELS cl dissipative pressure. 

Following [31] the equation for the evolution of the temperature of the 



radiation fluid is 



T 9H\ Op/dT 

T ~ Tdp/dT 6H dp/dT (ZZ) 
Inserting @, @ and (0) in (||) we obtain 

z[ _ 12K 2 (z) + 3zK 1 (z)-Bz 2 

z 12K 2 (z) + 5zK 1 (z)+z 2 K (z) ( ' 

where B = QHQ/Aq and ' = d/Hdt . From this equation it is easy to see 
that z' is negative for large z, which is in accord with the reheating scenario 
of above (see Fig. 2). 

We can explicitly obtain the time dependence of this temperature near 



t\ and t 2 . Expanding (18) about t = t\ and ( |2l|) for z — > oo, we obtain 
T oc t — t\. In the opposite limit (i.e. t — > t 2 and z — > 0) it follows that 
T oc (t 2 — t)~ l l 2 . Thus -as found in similar scenarios (3^]- the production 
of relativistic particles at the final stage of the black holes evaporation is 
accompanied by a huge increase of the temperature of the cosmic fluid. As 
for the temperature of the black hole component we must say that this is 
essentially zero since the black holes behave as a dust fluid (only that its 
"particles" emit radiation and (simultaneously) absorb the ambience fluid). 
Nevertheless, one may choose to ascribe a temperature to each individual 



mini-black hole of mass M by the Hawking relationship T^ t oc M -see [33]. 
The evolution of this temperature is governed by the sum of two terms. One 
of them comes from the black hole evaporation oc M~ 4 [33], and the other 



one comes from the accretion. The latter term is more complicated and not 
of much interest for our purposes here. In any case, the fate of the black 
holes (assuming they do not leave any stable relic behind) is their complete 
disappearance by yielding their whole mass to the radiation fluid. 
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The entropy production per unit volume in the latter is given by p(| 



CT 7o 4 T 

It has the limiting behavior 



vr 2 9-y 2 C 3 
S=^ = ~Vk (24) 



3 7 r(l -e~ x f 
m(3e- x -2) 4 (t 



- ^ „ ' , - , t-ti (25) 



where x = (£2 — ^i)/ T - Thus we see that the entropy production rate is 
very high at the beginning of the evaporation, while it decreases sharply 
at the final stage of this process. At first sight this may seem counter- 
intuitive if one has in mind that the final stage of black hole evaporation is 
explosive (when one adheres to the Hawking picture as we do). However, 
our assumption of £ = constant implies that the aforementioned accretion 
renders the evaporation rate much milder. The net rate of radiation particle 
production per mini-black hole and unit of volume roughly varies as (p+p)^ 1 , 
where in this case p and p refer to the radiation fluid only [34]. 



We note that the opposite process, i.e. the one in which the accretion of 
radiation by the black holes overpowers the evaporation of the latter, and as 
a consequence the radiation is entirely eaten up by the black holes, is ruled 



out by the second law of thermodynamics [32]. 



2.4 Variable 7 and r 

Assuming now that both 7 and r change in time, with a time scale for 7 
much larger than r, we may neglect the term r 2 7 in equation (|i~4|), whence 
we are left with D = — (1/4) ( 7 /t)' -ff 3 . Thus, the exponential inflation is 
stable provided 7/r is an increasing function, and whenever it begins to 
decrease, the de Sitter stage ends. 

The solution H = corresponds to the Minkowski solution. To analyse 
its stability we linearize equation @ about it, 

H + t^H - ^C,t~ 1 H = (27) 
The roots of the characteristic polynomial are 
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A± = ^ (-1 ± v/TTeCT) (28) 

Since £ and T are positive definite quantities it follows A_ < < A+ and 
therefore the Minkowski solution is unstable. There exists, however, a one- 
parameter family of solutions that approaches a flat spacetime solution at 
large times. There is in addition a one-parameter family of solutions that 
starts from a Minkowski spacetime in the far past and evolves towards a 
stable de Sitter solution. The time-reversal of the latter is nonsingular and 
corresponds to a spatially flat universe. 



3 General case 

In this section we consider the de Sitter solution (|l^) and assume that 7, r 
and £ are arbitrary functions. This more general situation may occur during 
the decay of massive particles into lighter ones and also during the decay 
of four-dimensional fundamental strings into massive and massless parti- 
cles -admittedly this second possibility is more speculative. The differential 
equation to solve is 

7 + r 7 = (29) 

We introduce a new dimensionless independent variable dr/ = di/r. The 
reason for using a dimensionless equation is that the equilibrium point of 
the differential equation (29), describing exponential inflationary models, 



will represent self-similar cosmological models. Then the general solution 
reads 
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^~"(/ dr/C^ + C) (30) 

where C is an arbitrary integration constant. We will use this result to 
present two simple models. First we consider that £ = Co> a constant. 
Thereby 

7(77) = ^(Co + Ce-") (31) 

A natural choice is r/ = z^lnt, v being a positive constant, that represents a 
linear dependence between r and the cosmological time t. Then we obtain 
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7 (t) = 70 1 + - (32) 



with 70 = Cq /Hi. This expression is monotonic increasing (decreasing) for 
C < (C > 0). For the case that C = (tih) v / (4% - Zt\) and 7o = 
(4*2 — 3ii)/(3(t2 — *i))) this model describes decay of massive particles into 
radiation beginning at ti and finishing at ti- 



Now inserting (32) into ( |15|) we obtain an expression for the dissipative 



contribution to the speed of sound 



This is a monotonic decreasing function that in the limit t — > oo behaves as 
v 2 ~ u/(3Hit). Thus we require that v(ti) < 1. 



Considering now that ( = Q^e v , the expression ( |30D gives 



7(^) = ^(^ + C) (34) 

In order to simplify the calculations we restrict to constant v; in this case 
we obtain 

&t = -^ln\ V + C\ (35) 
In order to have At > 0, we choose r\ + C > 0; then 

7(At) = 70 exp (C + 3#it> 2 At - e 3 ^^ A *) (36) 

This expression has a maximum at At = 0, and assuming that 7(0) = 4/3 
we find that there is a phase of decay into radiation starting when 7 = 1 at 
3H\v 2 At = —0.8678 that reaches a relativistic gas state at At = 0. After 
that there is a condensation phase back into nonrelativistic matter that ends 
at 3H\v 2 At = 0.6736 when 7 returns to 1. A scenario compatible with the 
latter phase is the quantum tunneling of radiation into black holes [35]. 
This may arise very naturally because of the instability of the hot radiation 
against spontaneous condensation (It is altogether different from the 

whole disappearance of the radiation by black hole accretion) . 
During this period both the viscosity coefficient 

C(t) = Co exp (C-e m ^ M ) (37) 
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and the relaxation time 



r(t) = -j^e"^* (38) 

are monotonic decreasing functions. 

Again we may interpret this behavior in terms of a two-fluid model, 
where the viscosity coefficient arises because of the particle production 
process from the decay of massive nonrelativistic particles into light ones. 
Shortly after the beginning of the decay the particle production rate is large 
and the energy density of the fluid becomes dominated by the light compo- 
nent. Later on, as the decay rate slows down the effect of adiabatic dilution 
by the fast exponential expansion of the universe turns out to be more im- 
portant. Therefore the nonrelativistic particles dominate again, since their 
number density goes down as a~ 3 , while the relativistic component goes 
down at the faster rate of a -4 . 

4 Conclusions 

With the help of the Lyapunov method we have studied the stability of 
cosmic inflationary expansions driven by a dissipative fluid whose transport 
equation is of causal type. We required a slow transition from the symmetric 
to the broken phase. This transition may be even quasistatic. Exponential 
inflation occurs during the transition whenever the ratio £/ (7 + rj) remains 
constant, or at most varies slowly in time. The behavior of the scale factor 
changes gently from exponential expansion to Friedmannian (a(t) oc t n , 
with < n < 1) once the heavy particles have either decayed or become 
sufficiently diluted, in this way rendereing the viscosity negligible. 

We have found that the de Sitter solution is asymptotically stable for 
a wide set of reasonable fluid quantities of the dissipative cosmological 
medium. Thus the condition for inflation appears natural, and the ex- 
tremely fast supercooling followed by an intense reheating, so frequent in 
the literature, is avoided. 

To obtain analytic results for the continuous decay processes of decay of 
mini-black holes we have modeled the fluid by a Boltzmann gas. Thus we are 
able to obtain the time evolution of the temperature from non-relativistic 
to the ultrarelativistic regime reflecting in an increase of the adiabatic index 
from 1 to 4/3. In this way a reheating phase is shown to occur simultaneously 
with exponential inflation. 
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Concerning the structural stability of the de Sitter solution relative to 
the spatial geometry, we have neglected the effect of curvature in compar- 
ison with the energy density and effective pressure during an exponential 
inflationary stage. In the case of non-causal bulk viscosity, where n is al- 
gebraically determined by H, it is relatively straightforward to investigate 
this question [37, ||, [|. By contrast, in the causal theory tt is no longer al- 
gebraically determined by H but satisfies a transport equation that couples 
it differentially to the expansion. Curvature introduces the scale factor ex- 
plicitly into the Friedmann equation, and makes it much harder to decouple 
the equations. Thus the effect of curvature in the causal case is far more 
difficult to determine in general. This will be the subject of future work. 
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Appendix 

The stability of a solution of equation (||) can also be studied by lineariza- 
tion about this solution. However, this procedure leads to a linear differen- 
tial equation for the perturbation with time-dependent coefficients, through 
three unspecified functions: CW> 7(0 an d 7"(i). Consequently it is extremely 
difficult to determine the behavior of the perturbation for large times, as it 
requires to calculate the perturbation as a functional of these coefficients. 
This is why we have chosen the second method of Lyapunov. This method 
has led us to some general qualitative results about the stability of the de 
Sitter solution. 

For a time-dependent Lyapunov function we have used the following 
theorems [plfl : 

1. If a function V exists which is defined and whose derivatives V is 
a semidefinite function whose sign is contrary of that of V , then the 
solution x = of 

x[ = fi(t,x lr -- ,x n ) , i = l,2,---,n, (Al) 
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is stable. 



2. If a function V exits which is definite, and has an infinitesimal upper 
bound, if the derivative V is also a definite function whose sign is 
contrary to that ofV, then the solution x = o/(Al) is asymptotically 
stable. 



which clearly satisfies the hypotesis of the theorem provided that (/t is 
bounded (see equation (15)), in a neighbourhood of the de Sitter solution 
Hi = C/(7 + T7). 
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Figure Captions 

Figure 1. Potential V{H) of the equivalent mechanical system defined by 
equation (|l0|). 

Figure 2. Plot of the dimensionless time derivative of the inverse tempera- 
ture z' for B = 0.1. 
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